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Abstract. We study the phase diagram of a microscopic model for the 
superconducting iron arsenides by means of a functional renormalization group. Our 
treatment establishes a connection between a strongly simplified two-patch model by 
Chubukov et al.[Tj and a five-band-analysis by Wang et al.f^. For a wide parameter 
range, the dominant pairing instability occurs in the extended s-wave channel. The 
results clearly show the relevance of pair scattering between electron and hole pockets. 
We also give arguments that the phase transition between the antiferromagnetic phase 
for the undoped system and the superconducting phase may be first order. 
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1. Introduction 

The discovery of high-temperature superconductivity in the iron arsenides of the 
llll[3lll] and the 122 [Sj famihes with critical temperatures up to = 56K has triggered 
enormous research efforts |6]. From the available experimental data and from density- 
functional theory[7J it became clear rather quickly that the superconducting pairing is 
most likely due to electron-electron interactions. Furthermore, most evidence points to a 
weaker local correlation strength than in the high- Tc cuprates[H]. Hence, weak coupling 
theories should actually work well in order to understand the phase diagram of these 
materials. 

Very early, first random-phase-approximation (RPA) treatments indicated that the 
predominant pairing instability should be in the extended s-wave or s±-channel with a 
pairing amplitude of opposite sign on electron- and hole-like Fermi surfaces[9]. However, 
a more complete theory should also account for vertex corrections and competing 
channels. Renormalization group (RG) methods go beyond RPA and sum all one-loop 
diagrams on equal footing up to infinite order. In two inspiring papers [2], the Berkeley 
group has already shown that the gross picture from the RPA is reproduced by the 
RG approach for two different descriptions [10] of the band structure, and supports the 
extended-s pairing state. However, their treatments of two- and five-band models also 
unveiled some differences between the predictions of the models, e.g. regarding the 
order parameter in the spin- density- wave (SDW) state. The complexity of the multi- 
band models indicated by these works, in particular due to the rich phase structure 
of the interactions in the Brillouin zone when the band language is used, definitely 
calls for other independent studies which show that the results do not depend critically 
on the details of the model and approximations used. Another, much simpler RG- 
like approach was recently presented by Chubukov and collaborators [1]. They used a 
strongly simplified electronic structure with one hole-pocket and one electron pocket, 
both of circular shape. Furthermore, they assumed that the interactions near the Fermi 
surfaces only depend on the pocket and not on the precise location of a wavevector on 
the respective pocket. Also in this model, for quite general repulsive interactions, the 
leading pairing instability turns out to be in the extended s-wave channel. 

Our study below aims to interpolate between the Berkley group and Chubukov et al. 
Like in the simplified treatment by Chubukov et al.[T], we define the bare interactions 
in fc-space, only dependent on whether the wavevector is near the F or near the Ap- 
point. In contrast to the work of the Berkeley group [2J, where the interactions are 
defined in real space, this will allow us to assess the effect of the wavevector-dependent 
matrix elements one gets from the transformation from real and orbital space into the 
band language. As the underlying band structure might be sensitive to details of the 
DFT scheme used to obtain it, this is an important check for the robustness of the 
results. It brings us closer to answering the question what the minimal model for the 
iron arsenide superconductivity is. Opposed to [1], in our functional RG additional 
wavevector-dependence around the Fermi surfaces can develop during the flows. This 
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could possibly lead to different pairing states, but turns out to be rather uneventful in 
our calculations. Furthermore, we use a more realistic band structure with two hole-like 
and two electron-like Fermi surfaces, such that we can study the effects of the nesting 
mismatch and doping in more detail than it was done in Ref. This way we can 
get estimates for the energy scales of SDW and SC ordering, and get out some more 
quahtative trends. 

This paper is organized as follows. In Sec. 2 we outline the functional RG approach 
as a tool for detecting instabilities in a Fermi liquid. In Sec. 3 we briefly summarize and 
comment on the two-patch model used by Chubukov. In Sec. 4 we present our extended 
modeling of the system, and in Sec. 5 we discuss the fRG results. Sec. 6 contains a 
summery and discussion. 

2. Functional RG 

The renormalization group is a well-known general concept to derive effective theories, 
e.g. at long length scales or for a low-energy window of a given many-particle system. For 
weakly coupled fermion systems, one is mainly interested in the effective interactions 
near the Fermi surface, as they contain information about possible symmetry-broken 
or non-Fermi-liquid ground states. Here, RG approaches to interacting fermions are 
more unbiased than diagrammatic summations in a particular channel, as competing 
channels are treated on equal footing. The self-energy near the Fermi surface is also of 
high importance, but is typically harder to obtain. Moreover, many trends in the self- 
energy, like gap-openings are "fore-shadowed" in the effective interactions. To compute 
the effective interactions near the Fermi surface, many recent works use the RG flow 
equations for the effective action or one-particle irreducible vertex functionsfTT l [T2l [T3] . 
that avoid some complications of other straightforward adaptions of Wilsonian RG 
for interacting fermions |14j. These RG schemes are commonly dubbed functional RG 
(fRG) as they aim at keeping as much of the wavevector- and possibly also frequency- 
dependence of the vertex functions, i.e. describe a flow of coupling functions rather than 
a flow of a finite number of coupling constants. 
We study a model described by an action, i.e. 



s = ±1/2. Q{k) denotes the quadratic part of the action, here given by Q{k) = 



+ Jn T^^'k'^i ^(k, k', k + q)i>k+q,si'k'-q,s'i'k',s'i^k,s ■ ^"^^ 

Here ip{k, s) and iplk, s) are Grassmann fields representing fermions with wavevector k, 
band index b, Matsubara frequency ko (we write k = {ko,k,b)) and spin projection 




iko,k 
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C^[e{k)] is very large for \e{k)\ < A and C\[e{k)] = 1 for |e(A;)| > A such that modes 
below A are not integrated over in the functional integral. In practice, one mainly needs 
the inverse C^^[e(fc)], which is chosen as a smoothened-out step function for numerical 
treatments. The Green's function as the inverse of Q\ is suppressed for modes with 
|e(/c)| < A. Here we neglect any self-energy corrections. For spin-rotational and U(l) 
invariance, the antisymmetric four-point interaction vertex for incoming particles fci, Si, 
^2, S2 and outgoing particles fcs, S3, ^4, S4 (the quantum numbers of particle 4 are dictated 
by the conservation laws, except for the band index 64 of the second outgoing particle) 
can be expressed in terms of a spin-independent coupling function V{ki, k2, k^, 64) which 
still carries the symmetries of the lattice [12]. The central RG differential equation 
governing the flow of the effective interactions V^i^ki, ^2, k^, 64) is|12] 

-^VAiki, ^2, ks, 64) = Tpp,A + T/^^A + Tp- (3) 

with the one-loop particle-particle contributions Tpp\ and the two different particle-hole 
channels A '^ph \^ where 

^P,a(^1) ^3; ^4) = 

- dkY^ VA{ki, k2, k, b') U^v{k, -k + ki + ^2) VA(fc, -k + ki + k2, fcs, b^) (4) 

b' 

'^ph,k{^Ii ^2; ^3, ki) = 

- dk'^ [■~2Va(A;i, k, k^, b') L^^h'{k, k + ki - k^) V^ik + ki - k^, k2, k, 64) 

b' 

+ Va(A;i, k,k + ki- k^, b') L^^b'ik, k + ki - k^) V/y{k + ki - k^, k2, k, 64) 
+ Va(A;i, k, k^, b') Lbfi'ik, k + ki - k^) Va(A;2, k + ki - k^, k, 64) (5) 

^H,a(^15 ^2; ^3; ^4) = 

- / dk'^ VAiki,k + k2- k3,k,b') Lb,b'ik,k + k2 - k3)VAik,k2,k3,b4) 

b' 

In these equations, J dk = Tj^iko I (C^' band index b' of the second internal 

line is to be summed over, while its wavevector and frequency are determined by 
conservation laws, e.g. as —k + ^1 + ^2 in Eq. HI The product of the two internal 
lines in the one-loop diagrams is 

Lb,b'{k,k') = SA{k)G^^\k') + G^^\k)SAik'), (7) 

with the band indices b, b' belonging to k and k' and the so-called single-scale propagator 

d 



(6) 



SA{k) = -G^l\k) 



Gfip). (8) 



Here, G^l\k) denotes the full Green's function at the RG scale A, but as stated above 
in this work we do not include self-energy corrections. Equations ([3]) can be integrated 
numerically for initial values corresponding to the bare interactions. At low A and low 
temperatures T, the flow will typically lead to strong coupling, i.e. a certain part of the 
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V\{ki, /c2, k^) will flow to values much larger than the bandwidth. Then the flow has to 
be stopped as the approximations (neglect of selfenergy, truncation of the hierarchy of 
the flow equations) will break down. Nevertheless, the scale where the flow to strong 
coupling occurs gives an upper bound for the energy or temperature scales of ground 
state ordering. Most importantly, the analysis of which channel diverges most severely 
indicates the type of ground state state ordering. 

The fRG equations stated above are valid at any temperature. Estimates for 
critical temperatures of long-range-ordered phases can be obtained by determining 
the temperature above which the flow remains finite down to zero scale. However, 
in standard cases, this value for Tc can (up to a factor of order one) also be inferred 
from the energy scale for the flow to strong coupling at T = 0. For instance, for a 
reduced Cooper problem with attractive interactions, we would get Tc = 1.14Ac where 
Ac = Wexp{—l/g) with the bandwidth W and the dimensionless coupling constant g. 
Hence, in the following we only compute Ac(T = 0). Furthermore, the approximate flow 
equations used here do not capture the destruction of long-range order in two dimensions 
by soft modes in the case of continuous symmetries. Yet, in the present study, this is 
not a problem, as in the "back of our minds" we actually want to study quasi-2D 
systems with some coupling in the third direction. There this destructive mechanism by 
Goldstone-modes is not effective, and spontaneous breaking of continuous symmetries 
is indeed possible. 

3. Two-patch model 

In this section, we first describe the most drastic but still meaningful simplification of 
the functional RG for the pnictide problem. Photo-emission experiments and density 
functional theory on the iron pnictides reveal (at least) 4 Fermi surfaces in the Brillouin 
zone[6j. Using the unit cell with two iron and two arsenic atoms, one has two nearly 
circular hole-like Fermi surfaces around the F = (0, 0)-point and two more deformed 
electron-like Fermi surfaces around the M = (vr, 7r)-point (for details see Fig. [1]). The 
simplified model by Chubukov et al. [1] reduces these Fermi surfaces to two circles, one 
hole-like Fermi circle around the F- and one electron-like circle around the M-point. In 
a further reduction step, the interactions between the quasiparticles states near these 
Fermi surfaces were taken to depend only on whether the excitation occurs near the hole- 
like or the electron-like Fermi surface. For the particle-hole-symmetric case this leads 
to 4 independent couplings gi, . . . ,g4 which are indicated in Fig. [2](b). In comparison 
with Ref. |T] we have interchanged gi and g2. The RG flow, upon integrating out the 
single-particles state with decreasing energy scale A, gives rise to the equations[l]: 



gi = 2gi{g2 - gi) (9) 

92 = gl + gl (10) 

^3 = -2(73(74 + 2(73(2(72- (71 ) (11) 

9^ = -9l-9l (12) 
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Here the notation gi = is used, and p(0) is the density of states at the Fermi level 

per Fermi circle and spin component. The general flow for repulsive bare interactions 
leads to strong couphng, i.e. the g^s diverge at a nonzero scale Ac. g2 and g^ flow to +00 
while goes to —00. gi diverges less strongly. If one now considers how the effective 
coupling constants for various symmetry breaking channels flow, one finds that several 
channels develop instabilities at the same critical scale. For example the SDW coupling 
constant g2 + g^ diverges, as does the coupling constant in the s±-pairing channel, i.e. 
g3 — g^, and some other combinations. Since for repulsive interactions, the coupling 
constant in the SDW channel is always larger than in the pairing channel, Chubukov 
et al. interpret this flow as indicative of a SDW instability. While we believe that this 
is physically correct, it should be noted that in one-dimensional systems like two-leg 
Hubbard ladders one finds similar multi-channel instabilities with a joint divergence of 
SDW and unconventional pairing channels. In these systems , however, the ground state 
does not develop SDW order but maintains a spin gap, and no power-law correlations 
occur. In physical terms, the competition of several ordering tendencies does not result 
in a clear "winner" with some kind of (quasi-) long-range order, but in a gapped spin- 
liquid-like ground state that embodies various types of strong short-range correlations. 

In the efforts to understand the one-band Hubbard model on the two-dimensional 
square lattice near half-filling, many authors have studied the so-called two-patch 
model[T5]. Here one assumes that the Fermi surface is near the van Hove points at (0, vr) 
and (tt, 0), and that only these fc-space regions are important for the breakdown of the 
Fermi liquid one is interested in. Again, in the standard treatment, the interactions 
are taken to depend only on the patch index, and not on the precise location of 
the wavevector inside the patch. This way one has to consider again four coupling 
constants gi to g^, defined very analogously to the pncitide toy model described above. 
It turns out that the RG equations for the one-band Hubbard model in the two-patch 
approximation are just the ones of the pnictide two-patch model described above, with 
the slight difference that now the density of states actually diverges at the Fermi level 
placed at the van Hove energy. However, on this level of approximation, the density 
of states only regulates the magnitude of the critical scale, and does not change the 
interplay of the various couplings and symmetry breaking channels. Again the flow is 
given by g2 and g^ 00, while g^ —00, and AF-SDW, (i-wave pairing and d-density 
wave susceptibilities diverge together. In the cuprates and in analogy with the ladder 
problem mentioned above, this picture was interpreted as an indication for a complex 
"mother state" with several competing channels embedded, with possible relevance 
to the pseudogap of the underdoped cuprates [Ml |17]. Note that very recently, the 
interplay of several ordering tendencies in the model and approximation by Chubukov 
et al. Ilj has been understood [HI [19] as a S0(6)-symmetric flxed point of the RG flow 
with degeneracy of pairing, SDW and orbital density wave channels. It would be very 
interesting to learn more about properties like quasiparticle spectra in the vicinity of 
this unstable flxed point. 

For the pnictides, there is at present no indication to argue for a complex spin-liquid 
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state. Hence one might ask if the minimal two-patch model overestimates the channel 
coupling and suggests a more sophisticated picture than what is actually found in the 
experimental system. In terms of the SO(6)-symmetry|18]. the flow seems to veer away 
from the fixed point relatively early and is instead attracted by other fixed points with 
only one type of ordering. 

For the pnictides, there is at present no indication to argue for a complex spin-liquid 
state. Hence one might ask if the minimal two-patch model overestimates the channel 
coupling and suggests a more sophisticated picture than what is actually found in the 
experimental system. Indeed, in the following, we will show that a generalization to 
more than just four coupling constants results, depending on the band filling, in two 
quite distinct regimes with only one predominant instability. Furthermore, one could 
ask if the simplification to the 4 coupling constants with no dependence is legitimate or 
if additional, possibly relevant wavevector-dependences arise in the fiow. Below we will 
see that is not the case. 

4. Model and FRG 

To model the band structure of the iron pnictides, we employ a four-band model in 
the folded Brillouin zone (BZ) proposed by Korshunov and Eremin [20j. The free 
Hamiltonian reads as 

H, = -J2^'n,Ak) -J^^ldlmM, (13) 

k,i,a k,i,cr 

with band indices i = ai,a2,/5i,/^2, onsite energies and an electronic dispersion 
given by 

^F'"^ = t"^'"^(cos + cos ky) + ^2^'"^ cos k^ cos ky, 

= t{''^' (cos k^ + COS ky) + 1'^^ COS ^ COS ^. 

Using the condensed notation (e*, t\, tg), we assume parameter values of (—0.60, 0.30, 0.24) 
and (-0.40,0.20,0.24) for the ai and 02 bands, respectively, and (1.70,1.14,0.74), 
(1.70, 1, 14, —0.64) for the [3i and [32 bands (all values are in units of eV). The doping 
is controlled by adding a chemical potential term to Eq. f|T3l) . A properly adjusted 
then yields two circular hole pockets around (0, 0) and two elliptic electron pockets near 
(tt, vr) in the folded Brillouin zone as depicted in Fig. [UJ^b). For the undoped situation, 
the areas of electron- and hole-pockets are the same, and the total electron number per 
spin projection and unit cell is = 2, or n = 4, if we sum over the spins. As we have 
2 Fe ions in the unit cell, a doping of x electrons per Fe for the 1111 compounds gives 
a band filling n = 4 + 2x. The band dispersion of the model Hamiltonian ( |T3l) . as well 
as its Fermi surface topology at half-filling (/i = 0), is shown in Fig. [H Here, it is 
clearly visible that hole and electron pockets are almost perfectly nested for Q = (tt, tt) 
which causes a logarithmic divergence in the particle-hole channel and, thus, drives the 
(tt, 7r)-SDW instabihty. 
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(0,0) (7r,7r) ^ (tt.O) (0,0) -tt tt 

Figure 1. (color online) (a) Energy dispersion along the main symmetry points of the 
first BZ. (b) Fermi surface topology at half-filling consisting of two hole pockets (blue) 
aroimd F = (0, 0) and two electron pockets (red) near M = (tt, tt) in the folded BZ. 

Next, we include electron-electron interactions consisting of intra- {g^) and 
interband interactions {gi and (72) as well as an interband pair hopping (g^) in direct 
analogy to the two-patch model by Chubukov et al. [1] discussed in the previous section. 
Summing up these contributions, the interaction part can be written as 

k,k',q 
<t,<t' 

+ f E diAk + q)4Ak' - q)d,Ak')daAk) + h.c. 

k,k\q 
<t,ct' 

+ f E ^-(^ + - ^d^Ak')dpM + h.c. (14) 

k,k' ,q 
<y,<yl 

+ f E ^-(^ + ^'^'o.A^' - ^d^Ak')daM 

k.k',q 
<t,<t' 

+ f E + - Q)d^Ak')dpAk), 

k,k' ,q 

where the sum over all band index combinations with a G {ai,ai} and [3 G {/3i,/32} is 
implicitly assumed. For the numerical handling of the RG flow in Eq. ([3]), we discretize 
the Brillouin zone into segments (patches) as shown in Fig. [2](b), in generalization to 
the fRG-treatments of the one-band Hubbard model [2n [T6] . Each segment is confined 
within two neighboring dotted lines running from (tt, tt) to one of the corner points 
(tt ± TT, TT ± tt). If we consider an arbitrary momentum ki and a freely chosen band 
index 6j, the point ki?(/cj,6i) is defined as the intersection of the angle bisector lying 
in the same segment as ki with the FS sheet corresponding to band index hi. This 
provides AN different discretization points pictured as colored markers in Fig. EJ^b). 
Points 1-24 run counter-clockwise around the first electron pocket, points 25-48 around 





Figure 2. (color online) (a) Discretization of the Brillouin zone with a 96-patch 
system, (b) The four relevant scattering processes near the FS with coupling constants 
gi (interband), 172 and 174 (intraband) as well as an interband pair hopping g^. For 
each process gi, the two arrows denote the scattering of quasiparticles with fixed spin 
projection s on one line and s' on the other line. 



the first hole pocket, points 49-72 around the second electron pocket and points 73- 
96 from the second hole pocket. Since the leading part of the flow is given by the 
coupling function on the Fermi surface and at zero frequencies [I2], we approximate the 
effective coupling Va(A;i, k2, k^, 64) by V\{ki, k2, k^, 64) = V\{ki, k2, k^, 64), where we used 
ki = (0, kp'(A;j, bi), h.j). This implies that each momentum ki^ ^ can be varied within its 
segment, while leaving the coupling function Va(/i;i, ^2, ^3, ^4) equal to VA(fci, A;2, A;3, 64). 
Due to this approximation, it is useful to decompose the fc-space integration in Eq. ([3]) 
into a radial part with a constant coupling in the integrand and a subsequent averaging 
over all segments. Starting from an energy scale Aq of the order of the bandwidth and 
VAg(/ci, /c2, ^3, 64) determined by the bare interactions ^fi, . . . , 5^4, the flow equation can 
now be integrated down to a critical scale Ac where certain parts of VA^(fci, A;2, ^3, 64) 
start to diverge. Similar to the RG flow of the effective coupling function, we also 
compute the flow of certain susceptibilities in order to have a clear criterion for the 
leading instability. Therefore, we consider the coupling of external fields $sc and ^sdw 
to fermionic bilinears written as 

Hsc =^scy^Kc{k)(d\ ^ , -d\ dK ) 

k,s 
a. a' 

k,s 
a, a' 

with band indices a, a' running over the four Fermi surfaces. The form factor is taken to 
be band-independent, but in the numerical implementation for the pairing case, a and 
b are restricted to be either both on electron-pockets or both on the hole-pockets. For 
the SDW channel, one particle at a hole-like Fermi surface and one in at an electron-like 
Fermi surface. 
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These couphngs are renormahzed by one-loop corrections that involve the scale- 
dependent interactions. In the following, we will only consider static external couplings. 
Then we can set the frequency flowing through these loop diagrams equal to zero. This 
amounts to Qo = for the SDW case and zero total incoming frequency for the pairing 
channel. Using the notation k = (0, k, a) with a' for the second external band index, and 
k' = {ik'^, k', b) and b' for the second internal band index, the one-loop renormalization 
of these couplings is obtained by 

-^KdUk) = - J dk KdUk'WAik, k' , k' + Q, a')Lb,b'ik' , k' + Q) 

^h,,{k) =- j dk'K,{k')VA{k,-k,k',b')U^,ik',-k'), 

with Lfe^b/(-,-) defined in Eq. (JTj) and the summation over all involved internal band 
indices implicitly assumed. The one-loop flow of the susceptibilities is calculated in a 
similar manner, i.e. 

d f 

'J/^Xsdw = - dk hsdn,{k)Lbb'{k, k + Q)hsdw{k + Q), 
-^Xsc = - j dk hsc{k)Lw{k, -k)hsc{,-k). 

As initial conditions for the couplings, we choose hsdw{k) = 1 for the spin coupling and 
hl^(k) = 1, hlf{k) = cos{kx) + cos{ky) or h'^^{k) = cos{kx) — cos{ky) for the different 
pairing channels. The initial values for the susceptibilities are zero. It becomes obvious 
that certain classes of coupling functions, as for instance V{k, —k, k') and V{k, k', k'+Q), 
drive the pairing or, respectively, the spin susceptibility and hence indicate a possible 
instability in this channel. 



5. Flows to strong coupling and phase diagram 

In the following, we calculate the RG-fiow of the coupling function V{ki, ^2, /cs, b^) and 
of the susceptibilities described in Sec. 4. Using a 96 point discretization scheme, 
we examine the results for several fillings in the electron and hole doping regime. As 
bare interaction parameters, let us first choose gi = 92 = 9^ = 9i = OAeV. The flow is 
stopped if one element of V{ki, k2, /cs, 64) exceeds 25eV. Starting with the half-filled case, 
the coupling function V"(/ci, ^2, ^3, 64) in Fig. [3](a) with the third momentum fixed to 
point 26 and the fourth band index set to 2, clearly shows a vertical diverging feature. It 
corresponds to couplings V{k, k', k' + Q) with a fixed momentum relation of Q = {it, it) 
between the second and third momentum. As discussed in Sec. 4, such couphngs 
with fixed wavevector-transfer drive the spin susceptibility with ordering momentum Q 
and, thus, indicate a SDW instability. Correspondingly, the SDW instability in Fig. 
^d) diverges steeply at the critical scale, leaving behind the extended s-wave pairing 
susceptibility. 

Note that the instability described here only signals that the magnetic moments in 
neighbored unit cells are oriented in opposite directions. In order to infer information 
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Figure 3. (color online) (a-c) The coupling functions VA(fci, fc2, /ca, 64), discretized at 
96 points according to Fig. HJa) with the first outgoing wave vector fixed at point 
26 and the band index of the fourth momentum set to 2 corresponding to the first 
hole pocket with points 25-48. (a), (b) and (c) are calculated for electron dopings of 
X = 0.0, X = 0.08 and x = 0.16, respectively. The colorbars indicate the values of the 
couplings and for k2 = 2, and fcs — ^2 = (tTiTt). (d-f) Flow of the s^-wave (solid) and 
(tt, 7r)-spin susceptibility (dashed). 

of the ordering within the unit, we would have to define a real space picture, i.e. 
to determine the orbital and sublattice site content of the bands of our dispersion. 
The dominance of the SDW-susceptibihty is still found at small nonzero doping, but 
the critical scale goes down gradually with doping and the pairing channel becomes a 
competitor, as seen in Fig. Mh) and (e). 

For electron doping of x = 0.16 in Fig. [3](c), the leading divergence has changed 
significantly. Here, we find diagonal features which signify an intense scattering 
V{k, —k, k' ,b4) of total-wavevector-zero Cooper pairs. Moreover, there is a clear sign 
reversal in the pair scattering at point k = 24, k = 48 and k = 72, which represents 
the change of k from a hole-like to an electron-like FS and vice versa. If we consider 
{cos{kx) + cos{ky)) as the simplest lattice harmonic generating this sign change, the pair 
scattering is properly fitted by V{k, —k, k', —k') = —{cos^k^) +cos{ky)) {cos{k'^) +cos{ky)) 
with the incoming pairing partners k,—k on the same band a and the outgoing pair on 
band b. As can be seen in Fig. [3] c) there is also a growth of the scattering of interband 
pairs, e.g. with the partners on the two different electron pockets, that is however 
subleading, as the total wavevector of such pairs adds up to a small nonzero value. This 
difference is due to our patching scheme. As a word of caution, we cannot exclude that 
the interaction of zero-total momentum interband pairs might become as strong as the 
one of intraband pairs. This issue will be scrutinized in further studies. 

An alternative analysis of the pairing symmetry can be obtained by optimizing the 
coupling strength functional [IQl which in our case reduces to an eigenvalue problem for 
the pairing matrix Vpair{k, k') = V{k, —k, k', 64). Here, the eigenvalues and eigenvectors 
of Vpair represent the pairing strength and the corresponding form factor. For the largest 
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Figure 4. (color online) (a) Flow of the d— wave (dashed), s-wave (blue solid) and 
s^-wave (black solid) pairing susceptibility at 16% electron doping, (b) The form 
factor fsc{k) (black) determined by the eigenvector corresponding to the largest 
eigenvalue in magnitude of the matrix Vf^^{k, —k,q). For comparison, we also plot 
0.1(cos(fcj;) + cos(fcj,)) (red) evaluated at the discretization points in Fig. [51 



eigenvalue in magnitude, we depicted the corresponding eigenvector in Fig. 111(b) and also 
calculated the flow of the pairing susceptibility for different form factors and intraband 
pairs in Fig. |4](a). 

These results clearly prove an extended s-wave pairing symmetry in accordance with 
the previous RG studies [2l H]. Note also, e.g. by looking at the plots of the low-scale 
interactions in Fig. [3l that we find very clean SDW at a; = or pairing instabilities for 
X = 0.16, respectively. In these plots, only a single feature, i.e. either the vertical lines 
plus the related weaker horizontal feature of an effective long-range SDW interaction, or 
the diagonal feature corresponding to the pairing instabilities, are visible. This shows 
that the multi-channel instability found in the two-patch analysis overestimates the 
channel coupling, and Chubukov et al. [T] were correct in their interpretation as a SDW 
instability. Furthermore, near the critical doping where the pairing susceptibility begins 
to grow more strongly than the SDW susceptibility, there is little overlap between SDW 
and pairing features in the effective interactions. This is unlike the situation in the 
one-band Hubbard model in the saddle point regime [16] and is well compatible with 
a simple change in the ground state from one type of order to the other without any 
exotic intermediate phases. 

From the pairing form factor in Fig. 111(b) we can infer that the pairing is equally 
strong on all Fermi surface sheets. This is a difference to conclusions drawn from 
photoemission experiments that showed a smaller gap on one of the hole pockets [22]. 
Our results indicate that this disparity may be due to a more detailed structure of 
the interactions near the Fermi surfaces, which is not included here. As another 
difference to microscopic studies based on more realistic five-band models [2], [23] is the 
mild variation of the pairing amplitude on the electron pockets. In the work by Graser 
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et al.|23j, the pairing even changes sign on the electron pockets, however in a way 
preserving the extended s-symmetry. Again we suspect that the matrix elements from 
the transformation from the full band structure to the band picture included in these 
works and not in ours is responsible for this difference. 

Pairing components different from the extended s-wave channel flow much more 
weakly. In particular, the d-wave component found in several studies [21 [23] is practically 
invisible in the effective couplings, and its susceptibility only shows a mild rise at low 
scales. The RG-fiow for the hole doped case is not depicted here since the findings are 
quite similar as in the electron doped regime. 

Next, we determine the critical scale Ac as a function of doping and examine which 
of the susceptibilities Xsdw and xlt diverges first. The results are summarized in the 
phase diagram of Fig. [H Near x = 0, the critical scale is highest, while it decays 
away from x = 0. On the tails, we have marked the regions where the extended s-wave 
pairing susceptibility diverges more strongly than the SDW component. Even though 
the absolute value of critical dopings and critical scales strongly depend on the values 
of the initial interactions, we find an interesting trend in the data: If we compare the 
critical scales for the pairing instability at doping ±x, i.e. the same absolute value of 
hole or electron doping, the critical scale is higher on the electron-doped side. In the 
experiments, although more factors might influences the situation as well, the higher 
Tc's are likewise measured on the electron-doped side. In our theory, this trend is not 
easily understood in terms of averaged quantities. For example, the density of states 
near the Fermi level is lower on electron-doped side, and also the bare spin susceptibility 
at (vr, vr) is stronger on the hole-doped side. Hence, we believe, that it is essential to 
resolve the wavevector-dependence of the interaction as done by the functional RG in 
order to see this trend. 

Next, we describe the dependence of the pairing instability on the initial 
interactions. Some data is shown in Fig [5b). One clearly sees that the critical scale 
dramatically depends on g^, with g^ = 0, where we do not find any pairing instability at 
reasonable scales. This fact can be understood already in the two-patch model and is 
also stated in Ref. [1]. Our results show that the refined wavevector-dependence does not 
alter the situation, and that no other significant pairing instabilities occur if the extended 
s-wave channel is taken out by g^ = 0. Hence a certain amount of pair scattering 
between electron- and hole-pockets is essential for a sizable Tc. Furthermore, we observe 
a repulsive interaction gi works against the pairing instability, but is otherwise not of 
central importance. 

Regarding the transition between SDW and extended s-wave pairing, we observe 
in the plots in Fig. [3] for the coupling function that the SDW channel has only little 
overlap with the pairing channel. Hence, it is reasonable to assume that its flow and also 
the effective potential for the SDW order parameter will be only weakly affected by the 
pairing correlations at low T. Then the SDW order may reasonably well be described 
by a renormalized mean-field theory. Meanfield studies of density-wave instabilities in 
imperfectly nested situations quite generally give a regime with first-order transitions. 
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doping X [eV] 

Figure 5. (color online) (a) Phase diagram for different initial interactions. The 
crosses (red) and squares (black) indicate the SDW and superconducting transition, 
respectively. For the upper line we used 51 = <?2 = 33 = 54 = OAeV as initial 
interaction parameters and for the lower one gi = gs — O.OSeV, 32 = 34 = OAeV. 
(b) Effective coupling Va(Ai, fe, 48, 4) near the critical energy scale Ac for initial 
interactions gi = gs = OeV, g2 — gi — OAeV and electron doping x — 0.16. 

where above a critical doping x the first-order transition hne Tt{x) extends from a 
tricritical point above the second order transition line Tc{x) and drops to zero at some 
higher critical doping XcjSl]. In our case, we have to keep in mind that the functional 
RG in this version is not sensitive with respect to first order- transitions |24]. But 
as the channel coupling is rather weak, regarding the effective potential for the SDW 
order parameter, we expect a similar situation as in the meanfield picture. Then, the 
SDW phase would actually extend to higher doping than is depicted in Fig. [5](a), with 
the shoulders given by first order SDW-to-metal transitions at temperatures above the 
pairing scale. As the energy difference between metal and superconductor is small, it is 
very likely that also the superconductivity-to-SDW transition at low T is of first order 
as well. It will be interesting to check this scenario with more sophisticated calculations. 
Very recently, the meanfield phase diagram for competing extended s-pairing and spin- 
density wave has been worked out in a two-pocket model for the pnictides|25j. confirming 
the basic picture of a first order transition, with the possibility of an intermediate phase 
with incommensurate SDW phase coexisting with pairing for large differences of the 
energy scales for pairing and SDW. In our RG study we refrain from distinguishing 
between commensurate and slightly incommensurate order as the reduced wavevector 
resolution would not permit definite statements. 

6. Summary and discussion 

Our functional RG study of the four-band model clearly identifies SDW order near 
zero doping and extended s-wave pairing with a sign change of the pair amplitude 
between electron and hole pockets at nonzero electron and hole doping as the two 
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leading instabilities. Hence, the basic picture for the instabilities by Chubukov et al.[T] 
remains valid, even if one allows for more realistic non-spherical Fermi surfaces and 
more possibilities for variations of the coupling function within the Brillouin zone. The 
additional variation of the dominant interactions around the Fermi surface turn out to 
be mild, or irrelevant in the RG sense. According to our generalization, this two-patch 
physics connects smoothly to the more ambitious studies by the Berkeley group [2] who 
take into account the ab-initio band structure and start with multi-orbital interactions 
in real space. Compared to our study, this changes details of the band structure and 
adds additional wavevector dependence to the initial interactions. Yet, the leading 
instabilities turn out to be of the same type. 

The main reason for the extended s-wave pairing is identified as the Qs-like 
scattering between electron and hole pockets. This interaction also drives the SDW 
tendencies, as was emphasized recently by the Berkeley group [26]. This important role 
of the particle-hole channel with momentum transfer (vr, vr) is very similar to the weak 
coupling picture of the d-wave superconductivity in the one-band Hubbard model near 
half filling, which is possibly also relevant for the high-T^ cuprates. In this sense, and if 
our modelling indeed captures the relevant physics, the pairing mechanism in cuprates 
and iron arsenides can be considered to be of the same type. The difference in the pairing 
symmetry is due to the different location of the Fermi surface parts responsible for the 
pairing instability. Furthermore, in the one-band Hubbard model on the square lattice, 
the competition and coupling between the various channels can be more intense p^. 

Our study allows to obtain estimates for the instability scale as a function of doping 
and the critical doping where the SDW phase gives way to the pairing instabihty. Here, 
we find that the electron-doped side with smaller density of states near the Fermi surface 
shows the higher critical scale for pairing if we compare the same absolute value of the 
doping. This trend is also found experimentally. It shows that in anisotropic systems 
with electronically mediated pairing not only the averaged density of states decides 
on the energy scale for pairing. This scale can be a rather complicated function of the 
dispersion and interaction parameters. Notably, the observed trend cannot be explained 
by a higher degree of nesting on the electron-doped side. At least, if we take the bare 
particle-hole susceptibility at wavevector (7r,7r) as a measure for the nesting, it turns 
out to be smaller on this side with the higher pairing scale. 

Concerning the transition between SDW and extended s-wave pairing, the 
functional RG exhibits a continuous change of the instability with doping, with a 
smooth variation of the critical scale. As we do not include self-energy corrections, 
this scale should be considered as an upper bound for the breakdown of the metallic 
state, and a second order transition between SDW and extended s-wave pairing cannot 
be excluded rigorously. However, as the couping between the channels is rather weak, 
and as commensurate density-wave phases in imperfectly nested situation often end in 
first-order transitions, we have argued that a first-order transition between the SDW 
and the superconducting state is quite plausible. 

Finally, we can compare our results with relatively simple bare interactions defined 
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in the band picture with other approaches where the interactions are defined in the 
orbital real-space picture. We observe two differences, which indicate the influence of 
the 'matrix elements' set to unity in our study. One aspect is the question of competing 
pairing channels. The fRG with matrix elements [2] and also RPA-type calculations |23] 
flnd a c?-wave component that is only slightly weaker than the extended s-wave channel. 
In our calculations, the c?-wave component is much weaker, and does not appear as a 
competing channel. Further, photoemission [22] and other theoretical studies [23] obtain 
a different magnitude of the pairing gap on the two hole pockets. Our results do not 
indicate strong differences between the pairing strengths of the 4 Fermi surfaces. Again, 
it is possible that this effect is due to the matrix elements that give additional structure 
to the bare interactions. 
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